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ABSTRACT 

We determine the ensemble properties of z ~ 5 Lyman break galaxies (LBGs) se- 
lected as V-band dropouts to iab < 26.3 in the Chandra Deep Field South using their 
rest- frame UV-to- visible spectral energy distributions. By matching the selection and 
performing the same analysis that has been used for z ~ 3 samples, we show clear differ- 
ences in the ensemble properties of two samples of LBGs which are separated by 1 Gyr 
in lookback time. We find that z ~ 5 LBGs are typically much younger (<100 Myr) and 
have lower stellar masses (^10 9 Mq) than their z ~ 3 counterparts (which are typically 
~fewxl0 10 Mq and ^320 Myr old). The difference in mass is significant even when 
considering the presence of an older, underlying population in both samples. Such 
young and moderately massive systems dominate the luminous z ~ 5 LBG population 
(> 70 per cent), whereas they comprise <30 per cent of LBG samples at z~3. This 
result, which we demonstrate is robust under all reasonable modelling assumptions, 
shows a clear change in the properties of the luminous LBGs between z ~ 5 and z ~ 3. 
These young and moderately massive z ~ 5 LBGs appear to be experiencing their first 
(few) generations of large-scale star formation and are accumulating their first sig- 
nificant stellar mass. Their dominance in luminous LBG samples suggests that z ~ 5 
witnesses a period of wide-spread, recent galaxy formation. As such, z ~ 5 LBGs are 
the likely progenitors of the spheroidal components of present-day massive galaxies. 
This is supported by their high stellar mass surface densities, and is consistent with 
their core phase-space densities, as well as the ages of stars in the bulge of our Galaxy 
and other massive systems. With implied formation redshifts of z ~ 6 - 7, these lumi- 
nous z ~ 5 LBGs could have only contributed to the UV photon budget at the end 
of reionisation. However, their high star formation rates per unit area suggest these 
systems host outflows or winds that enrich the intra- and inter-galactic media with 
metals, as has been established for z ~ 3 LBGs. Their estimated young ages are con- 
sistent with inefficient metal-mixing on galaxy- wide scales. Therefore these galaxies 
may contain a s ignificant fraction of meta l-free stars as has been previously proposed 
for z ~ 3 LBGs (jJimenez fc Haimanll200fih . 
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ies: starburst 
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1 INTRODUCTION 

One of the fundamental open questions in cosmology is when 
did galaxies form their first generations of stars? Identify- 
ing and studying such galaxies are key steps towards un- 
derstanding the physical processes that drive galaxy forma- 
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tion. Probing the formation and early growth of systems 
similar to the Milky Way requires observations of galaxies 
when the Universe is still youn g. Specifically, theoretical cal- 
culations l|Mo fc White 2002) predict that the most rapid 
growth of galaxies with masses comparable to that of the 
Milky Way occurs at redshifts of approximately z ~ 5. This is 
not long after the end of the (complex) reionisation process, 
which recent t hree- year WMAP r esults indicate was under- 
way at z ~ 11 (| Alvarez et al.ll2006|) and was largely complete 
by redshifts z ~ 6-6.5 jFan et al.1 l200ll; iBecker et all [200ll : 
iMalhotra fc Rhoadsll2004 iFan et al.ll2005l ). Estimates of the 



ionising photon density in the early Universe suggests that 
the UV emission from currently known high redshift galax- 
ies, whether star-formation or AGN dominated, is insuffi- 
cient to have caused reionisation. However, these galaxies 
must have had an impact on the intergalactic medium at 
high redshift. Only through comprehensive studies of the 
physical properties (including masses, star formation rates 
and histories, and clustering strength) of galaxies that were 
in place at this epoch, can we accurately assess their con- 
tribution to the mass growth of galaxies like our own, and 
their effect on the gaseous intergalactic medium (IGM) at 
the end of reionisation. 

As part of an ongoing study of high redshift galax- 
ies, we have investigated the rest-frame UV-optical prop- 
erties of candidate LBGs at 4.6 < z < 6 selected as V- 
dropo uts using the now standard Lyman-break tech - 
nique jSteidel fc Hamilton! Il993l ; ISteidel et al] 1 19951 . Il999h . 
Through spectroscopic confirmation, we have successfully 
demonstrated the efficacy of this technique for unam- 
biguously identifying z ~ 5 galaxies from deep imaging 
surveys in the res t -fram e UV using 8m-class telescopes 
l|Lehnert fc Bremerl 120031 ). From these data it has been 
possi ble to determine the comoving density l|Bunker et al.l 
120041). the likely contrib ution to the end of reionisation 
IlLehnert fc Bremer|[2003l ). and the fraction of sources which 
host super-massive black holes l|Bremer et al.|[2004l ) of z ~ 5 
LBGs. However, the rest-frame UV data alone are thus far 
insufficient to accurately constrain the ages, dust content, 
star formation rates and masses of the LBG population at 
high redshift. 

Rather, accurately constraining these parameters relies 
on well measured rest-frame UV-to-optical SEDs. At z ~ 5, 
the rest-frame UV emission from ongoing star-formation is 
redshifted into the observed visible, while emission in the 
rest-frame visible to near-infrared from older stars, diagnos- 
tic of longer or earlier periods of star formation, shifts into 
the mid-infrared. Several fields now have excellent multi- 
wavelength data sets from both the Hubble and Spitzer 
Space Telescopes (HST and Spitzer, respectively), supple- 
mented by ground-based data sets, forming an ideal basis 
for selecting and studying samples of distant galaxies across 
the full rest-frame UV to visible wavelength range. 

The most detailed multiwavelength studies to date 
of 5 < z < 7 LBGs have been centered upon a small frac- 
tion of dropouts that benefit from amplification due 
to lensing and/or those detected with Spitzer-IRAC 



the ages of these systems. A strong detection with IRAC 
normally confirms the presence of a Balmer break or an 
intrinsically more UV luminous system. The derived prop- 
erties of these high redshift LBGs are fascinating; massive 
(few-several x 10 10 Mq) and strongly star forming systems 
with ages comparable to the age of the Universe at that 
epoch. The presence of such massive and evolved systems 
at high redshift, in place less than a billion years after the 
Big Bang, challenges the expectations from bottom-up hi- 
erarchical structure formation scenarios (however, also see 
iMcLure et al1l2006l ). 

These derived masses and ages of z > 5 LBGs are sim- 
ilar to the average properties of LBGs at redshift 3. Fol- 
lowing the pioneering work of Steidel et al., the last decade 
has seen intensive observational and theoretical studies on 
the properties of LBGs at z ~ 3, providing a wealth of in- 
formation about this abundant population of UV-bright 
galaxies at this epoch. These star-forma tion dominated sys- 
tems are seen to host strong outflows (|Pettini et al 1 l200ll : 
lAdelberger et alj|2003l ; IShaplev et"ai1l2003l) , t ypically have 
sub-s o lar-to-solar m e tallic ities (~0.3-lZ©; IPettini et al.l 
l200ll : IShaplev et all 120031 ) and possibly contain a sig- 
nificant component of m etal-free (population HI) stars 
ij Jimenez fc Haimanl 1200?} ). LBGs are shown to be highly 
clustered on both large and small s cales, the latter being 
indicative of common halo objects rtGiavalisco et al.l 19981; 
Adelberger et al.lll998l : lOuchi et al.1l2004bl . 120051 : iLee et all 
20061 ). Recent follow-up of z~3 LBGs with Spitzer has re- 
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the Balmer-break lies between the K s and the IRAC pass- 
bands, the IRAC data are highly effective in constraining 



vealed a division in the population between LBGs with 
significant dust-attenu ated star-forming r egions as well as 
the UV-emitting ones (iHuang et al 1l2005l). This division is 
also seen for LBGs at z ~ 1 l|Burgarella et al.ll2006l) with 40 
per cent being infrared bright. Recently iReddv et al.l ((2005n 
have reported on the pot ential relatio n betw een UV- and 
IR selected populations. iHuang et al.1 (|2005l ) suggest IR- 
luminous LBGs are the missing link between LBGs and 
sub-mm galaxies and are the progenitors of present-day gi- 
ant e l lipticals (see also | Adelberger fc St eidel 2000; Shu et all 
l200ll : lQuchi et al.ll2004bl : iRigopoulou et aLlboOrj ). 

But how do these properties compare to those samples 
of LBGs at higher redshifts? Do we see an evolution in the 
properties of similar LBGs at these two epochs? Because of 
the small number of sources investigated, and the inhomo- 
geneity of the selection criteria used, it has not yet been 
possible to consistently compare the derived properties of 
z > 5 LBGs directly with thos e of Lyman- b reak galaxies at 
z~3. In a prelimi nary study, Ando et alj l|2004l ) find that 
unlike z ~ 3 LBGs (jShaplev et al .1120031 ) the majority of lu- 
minous z ~ 5 - 6 LBGs have weak or absent Lyman alpha 
(Lya) emission and strong low-ionisation absorption lines, 
already indicating key differences between the populations 
at redshift 5 and redshift 3. We complement this compari- 
son by analysing the properties of a large and reliable sample 
of the most luminous LBGs (L > L*) at z ~ 5 which match 
the luminosity to which LBGs at redshift 3 have been spec- 
troscopically confirmed. This matched selection permits a 
direct comparison between the properties of similarly lumi- 
nous LBGs in place in a ~1.2 Gyr old Universe and those at 
redshift 3, ~1 Gyr later. This difference is longer than the 
typical duration of the UV luminous phase of a z ~ 3 LBG, 
thus we are not comparing the same galaxies, but galaxies 
at two epochs with similar UV-emission characteristics. 
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To address these issues, we present the results from an 
analysis of the properties of a sample of luminous V-band 
dropouts selected from the multiwav elength datasets of th e 
Chandra Deep Field South (CDFS. Iclacconi et al.ll200ll ). 
This field was chosen as the southern field of the Great Ob- 
servat ories Origins Deep Survey (GOODS, iDickinson et al.l 
2003a), for which public imaging data are available over 
ten bands covering the optical (HST-ACS), near- infrared 
(VLT-ISAAC) and mid-infrared (Spitzer-IRAC) wavelength 
ranges. In section [5] we describe our methodology and selec- 
tion criteria that define this robust sample of z ~ 5 LBGs. 
The stellar evolutionary modelling is discussed in section [3] 
The resultant physical properties of mass, ages, star forma- 
tion rate and extinction of the galaxies in this sample are 
detailed in Section and the comparison to the proper- 
ties of z ~ 3 galaxies is presented in Section Unlike sim- 
ilarly luminous z ~ 3 LBGs, and other z > 5 LBG samples, 
we find our sample to be dominated by young (<100Myr) 
and moderately massive (~ 1O 9 M0) galaxies. We discuss 
the implications of our findings in Section [6] 

We adopt the following flat cosmology throughout this 
paper: 7io=70kms~ 1 Mpc~ 1 , O m =0.3, riA=0. 7. All magni- 
tudes are based on the AB magnitude system (|Oke fc Gunnl 
1 19831 ). 



2 DATA AND SAMPLE SELECTION 
2.1 Selection Criteria 

We applied a Lyman-break colour selection to the pub- 
licly available HST/Advanced Camera for Surveys (ACS) 
imaging datasets of the CDFS obtained through the fil- 
ters F43JW_J J B} J _F606W_( F), F775W (i) and F850LP 
(z) |Ciavalisco et al.l l2004bl ). Specifically, we selected ob- 
jects with Vab-iab > 1.7 mag which were not detected 
(signal-to-noise ratio < 3) in the F435W (B-band) image 
(i.e. short-ward of the 912A Lyman-break at z ~ 5) ensuring 
that we only selected sources with a clear Lyman discon- 
tinuity in their emission. We then required that the ob- 
jects were reliably detected (signal-to-noise ratio > 5 in 
the a-band) in order to accurately measure the sizes of the 
UV emitting regions. Furthermore, as mentioned above, we 
matched our sample to the magnitude limit to which LBGs 
at z ~ 3 have been spectroscopically confirmed. This limit of 
^B~25.5 probes z ~3 LBGs brighter than ~ 0.4L* Z=3 
jSteidel et all 120031 '). At redshift z~5 ^ B ~25.5 corre- 
sponds to iAs~26.3mag, the magnitude to which we select 
our sample of z ~ 5 LBGs. This colour selection is similar 
to those previously used to se lect V^-band dropouts (or the 
analogo us J?- band dropouts i n iLehnert fc Bremerll2003l . but 
see also IVanzella et al.ll2006l : IStark et all 20061 for a discus- 
sion of alternative colour cuts) , which have been successfully 
spectroscopically confi r med to lie at z ~ 5 ( Giavalisco et al.l 



l2004al ; lYan et alJlioolil ; IVanzella et ai1l2006l ) 



2.2 Input data 

For each object we combined the ACS data with publicly 
available near- and mid-infrared datasets to construct multi- 
wavelength SEDs covering 0.45 to 8/^m (Figure [TJ. The 
ground-based J- and A" 3 -band data were taken with IS AAC 
at the Very Large Telescope (VLT) (jOlsen et al.ll2006l, and 



http:/ /www. eso.org/science/goods/imaging/products. html ) 



1 We explore how different assumptions in the modelling process 
affects the results of our analysis in Appendix A and the potential 
contribution from an underlying old population in Appendix B. 



and the mid-infrared data with Spitzer/IRAC at 3.6, 
4.5, 5.6 and 8.0/xm (Dickinson et al., in prep, and 
http : // dat a.spitzer . caltech . cdu/ popular/ goo ds/ , see als o 
brief descriptions in lYan et al.ll2004l and lStark et al.l 2006'). 

While we used the publicly-available reduced images for 
the ACS and ISAAC data, we performed post-pipeline pro- 
cessing on the GOODS-S/Spitzer data (from epochs 1 and 2) 
creating deep 1' x 1' mosaics centered on each high redshift 
candidate using the Spitzer/MOPEX package. The resultant 
drizzled mosaics have had instrumental effects and cosmic 
ray events removed. The mosaics were registered to the ACS 
astrometric reference frame using common detections to an 
accuracy of ~0.1". We verified our IRAC mosaic generation 
and photometric calibration by performing the same proce- 
dure to e xtract photometry for A and F stars identified in 
the field l|Groenewegen et al.ll20o3 ). Such stars rarely show 
an infrared excess or have strong absorption features in the 
mid-infrared, thereby enabling accurate predictions of their 
mid-infrared emission. The fluxes predicted by black-body 
fits to the templates of their spectral types were reproduced 
to within 0.1, 0.1, 0.2 and 0.3 mag of our aperture photome- 
try at 3.6, 4.5, 5.8 and 8/im, respectively, consistent with our 
adopted accuracies (see later for a description of our adopted 
photometric uncertainties). Additionally, we confirmed that 
our measured fluxes were consistent with those measured 
for our LBG candidates and the stars from the epoch 1 and 
2 GOODS Enhanced Legacy products (Dickinson et al., in 
prep, and http://data.spitzer.caltech.edu/popular/goods/). 

The deep GOODS/IRAC data have the largest spatial 
pixels and PSFs amongst our datasets and, as a result, in 
crowded regions the emission profile from a given source can 
overlap with those of adjacent sources. For those sources 
selected according to the above criteria, but blended with 
neighbours in the Spitzer data, we treated the mid-infrared 
IRAC fluxes as upper limits in the subsequent SED analysis. 

For each source that satisfied our Lyman break crite- 
ria, we built a 10-band SED with fixed diameter aperture 
photometry performed on each image (see Figure [TJ. Aper- 
tures, appropriate for the characteristics of each image, were 
determined based on enclosing as large a fraction of the 
source flux as possible while minimising the impact of the 
sky background noise and flux contributions from neighbour- 
ing sources. We used apertures with diameters of 1", 2", 
and 4.5" for the measurements on the ACS, ISAAC, and 
IRAC images, respectively. Small aperture correction fac- 
tors were applied to account for the fraction of the total 
flux not enclosed by our fixed-diameter apertures: 1.07 for 
the four ACS bands, 1.10 for the two ISAAC bands, 1.18 for 
the IRAC 3.6 and 4.5 fim bands, and 1.35 and 1.47 for the 
5.8 and 8 /jm bands. We determined the aperture correc- 
tions based on the point spread function (PSF) of each mo- 
saic, constructed from isolated, bright but unsaturated stars 
throughout the CDF-S and taking reference total apertures 
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Figure 1. (i) This figure shows 5" x5" images of a typical galaxy from our sample in the 10 filter-bands we consider. Six spectral energy 
distributions (SEDs) are shown in panels (ii-vii). The best-fit model is plotted over the data points. Panels (ii-vi) show SEDs for five 
galaxies from our robust sample of z~5 LBGs. The SEDs are arranged (from left to right then down) in approximate best-fit age order 
showing the emergence of the Balmer-break between the K s and 3.6/xm data points. The SED of a low redshift interloper is shown in 
panel vii for comparison. Panels iii, iv and vi are spectroscopically confirmed z ~ 5 LBGs. 



of diameters 6" for the optical and near-infrared data, and 
12.2" for the mid-infrared data. 

The formal uncertainty assigned to each photometric 
data point was derived from measuring the background 
fluctuations in the parent image using ~10 3 -10 4 synthetic 
apertures with fixed diameters laid down at random on 
regions that are empty of flu x from discrete sources (e.g. 
iForster Schreiber et al.ll2006al ). thereby accurately charac- 
terising the photometric uncertainty directly from the data, 
rather than calculating it from the pixel-to-pixel rms assum- 
ing uncorrelated Gaussian noise statistics. The simulations 
were performed for a range of aperture diameters. As a func- 
tion of linear aperture size N, the empirically derived rela- 
tionship shows that the photometric uncertainties lie above 
and grow faster with size than the a oc N dependence for 
pure Gaussian noise. However, for a given aperture size, the 
histogram of background fluxes is very well described by 
a Gaussian distribution (of which the dispersion is taken 
as the la uncertainty for photometric measurements in the 
given aperture size). A conservative absolute calibration un- 
certainty of ten per cent was additionally included in the 
final adopted uncertainty to comfortably account for the rel- 
ative uncertainties in photometric calibration across the ten 
bands. In our SED modelling, this additional uncertainty 
prevents the fits being driven by a few photometric points 



with very small errors. We note that adopting a less conser- 
vative value for this uncertainty would mainly act to reduce 
the range of acceptable older ages, which is most influenced 
by the IRAC photometry, and therefore would not alter our 
findings of young ages for the majority of our sources. 

2.3 Defining a robust sample of z~5 
Lyman-break galaxies 

Our initial rest-frame UV selection yielded a sample of 109 
z ~ 5 LBG candidates. This sample will include a fraction of 
spurious sources that are not LBGs (low redshift galaxies, 
quasars and stars) and also some sources with insufficient 
photometric constraints (regardless of their nature) to re- 
liably determine their properties. We have therefore culled 
the inital sample as follows. 

Assisted by SEDs extending to the observed-frame mid- 
infrared, we have reliably excluded 22 candidates from our 
sample which are low redshift galaxies (or interlopers). Syn- 
thesis modelling of the colour evolution of stellar popu- 
lations predicts that galaxies with z < 4 can also satisfy 
our UV selection criteria. These galaxies have intrinsically 
redder near and mid-infrared SEDs than true high red- 
shift Lyman break galaxies. The SED of an example in- 
terloper is shown in panel (vii) of Fig. [T] and is easily dis- 
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tinguished from the SEDs of high redshift LBGs (panel 
ii-vi). Furthermore, Figure [T] demonstrates how the addi- 
tion of the IRAC data has greatly enhanced our ability to 
screen for such low redshift galaxies which are generally 
the brightest sources in the IRAC wavelengths that sat- 
isfy our optical selection criteria. These interlopers mostly 
satisfy the criterion for being extremel y red objects (ERO, 
i A B -K a ,AB >2.48, iRoche et all 1200311 and IRAC-se lected 
EROs (iERO, m 3 . 6 »m,AB - zab >20. lYan et~ai"1l2004h . con- 
sistent with being galaxies at z~ 1-2. The importance of 
this screening is clear when one considers that most samples 
of systems satisfying the high-redshift LBG selection crite- 
ria distinctly lack confirmatory spectroscopy and several of 
the recently reported IRAC-detected modelle d z > 5 LBGs 
do not have confirmed spectroscopic redshif ts l)Egami et al.l 
120051 ; iMobasher et al.ll2005l; lYan et al.ll2005l) . In the case of 
the z ~ 6 LBG in IMobasher et al.1 (120051)7 an interloper solu- 
tion is more plausible iDunlop et all d2006h . 

Our initial optical selection potentially includes low 
mass stars within the Galaxy and QSOs. These are spatially 
unresolved in the ACS data are therefore straightforwardly 
excluded. We identified 19 such objects. These objects gener- 
ally have the brightest visible emission amon gst the systems 
that s atisfy our selection criteria (zab ~23). lVanzella et al.l 
(2006) obtained spectra of four of these objects with suffi- 
cient quality for them to be classified as stars, supporting 
the exclusion of unresolved objects from our robust sam- 
ple. Had we included them and subjected them to synthe- 
sis modelling, their best-fit models would imply implausibly 
high star formation rates and young ages. Thus, without 
screening for these objects, our derived ensemble properties 
would be heavily biased towards young ages and high star 
formation rates. 

Finally, as accurate estimates of the ages, masses and 
star-formation histories of our galaxies are dependent upon 
well-sampled SEDs, it is imperative to analyse the objects 
with the most robust photometry. Therefore, we additionally 
excluded the results from 47 sources with secure detections 
in only two or three bands (including i and z), with the 
remaining constraints being upper limits. There are insuffi- 
cient data to constrain the models to accurately constrain 
their physical properties, and these systems are flagged as 
having 'insufficient photometry'. A significant fraction of 
these systems are blended with a nearby source in the IRAC 
bands. Because blending is purely a random alignment of 
foreground sources with our high redshift candidates, there 
is no reason these should be different to the isolated candi- 
dates. Accordingly, we do not find any evidence for a sys- 
tematic difference in the derived properties of the blended 
and isolated sub-samples. We note that this step may pref- 
erentially exclude LBGs with high star formation rates, low 
masses and young ages as these are less likely to be detected 
in the IRAC bands. Massive, old LBGs selected at the same 
i— band magnitude as younger, lower mass systems are in- 
trinsically brighter in the observed IR and so require more 
of the photometric bands to be compromised in order to 
exclude them in this way. 

This culling process renders a final robust sample of 21 
galaxies which is the focus of our study. Since each step in 
the culling process (apart from the last) is unbiased with re- 
spect to the physical properties of our high redshift galaxies, 
the final robust sample is expected to be representative of 
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Figure 2. Histogram showing the 2-band magnitude distribu- 
tions of the full set of 109 objects satisfying the initial selection 
criteria, with histograms of the comprising sub-samples overlaid. 
The black/shaded histogram represents the final 'robust' sample 
of reliable LBG candidates. 



the bright, z ~ 5 galaxy population selected by the Lyman- 
break technique, with the possibility that the lower mass 
and younger systems may be underrepresented. The use of 
the longer wavelength near- and mid-infrared data have en- 
abled us to construct a uniformly selected rest-frame UV 
sample that is reliable. Figure [5] shows the magnitude dis- 
tribution of the robust sample in comparison to all objects 
satisfying our selection criteria and the culled members. As 
expected, we find that the Star/QSO candidates and low- 
redshift galaxies are in the most part the brightest sources 
that satisfy our initial selection criteria. The histograms of 
the robust sample and the LBG candidates flagged as having 
insufficient-photometry suggest that the former comprises 
the brighter members of the z ~ 5 population and the lat- 
ter the fainter members. Indeed 6 LBGs from the robust 
sample, and 9 candidates with insufficient photometry, have 
been spe ctroscopically confirm ed to be z ~ 5 LBGs (see Sec- 
tion E31 |WnieiiaeraD[2Q0i). 



EVOLUTIONARY STELLAR SYNTHESIS 
MODELLING 



Using a library of synthetic spectra |Bruzual fc Charlotl 
2003), we modelled the SED of each LBG candidate simulta- 
neously deriving photometric redshifts and the key proper- 
ties of age, extinction, star formation rate and stellar mass. 
We explored a range of star formation histories (constant 
star formation rate, instantaneous burst, and exponentially 
decaying star formation rates with e-folding timescales rang- 
ing from 10 Myr to 1 Gyr). The UV emission from galaxies is 
dominated by the population of short-lived massive OB and 
A stars. Models of continuous star formation enable the pro- 
duction of these massive stars even at old ages. In contrast, 
massive stars die away rapidly in an instantaneous burst 
model and conspicuous UV emission is only produced at very 
young ages. For a given initial mass function (IMF), model 
fits with these two star formation histories will bracket the 
possible age range of the galaxies. Our emphasis on the con- 
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stant star formation model therefore provides upper limits 
to the best-fit stellar ages. For this reason we concentrate our 
analysis on the results obtained for a constant star forma- 
tion rate. The effects of adopting alternative star formation 
histories on our results are discussed in Appendix lAl 

For all of our models, we used a Salpeter stellar IMF 
between 0.1 and 100 Mg. The actual IMF is unconstrained 
for our objects. A steep rise down to the lower mass cutoff 
likely is unrealistic in view o f the turnover below 1 Mr in - 
ferred for the local IMF (e.g. lKroupall200ll ; IChabrierll2003h . 
and recent comparisons of dynamical masses and photo- 
metric stellar masses suggest that Kroupa/Chabrier-type 
IMF s would be more appropriate at higher redshift as well 
(e.g. IForster Schreiber et all l2006bh . The main impact of 
changing the IMF on our SED modelling is on the derived 
masses and star formation rates which, for a given rest- 
frame V-band lu minosity . would be about a f actor of 1.4- 
2 lower with the iKroupal ll200lfl or IChabried <|2003h IMFs 
(see [B ruzual & Chariot 2003; see a l so the discussion by 
IPapovich et al.l l200ll . in the context of SED modelling of 
z ~ 3 LBGs). At z > 2, and more so at z ~ 5, top-heavy IMFs 
with a larger proportion of massive (> 10 Mg) stars, or 
even extremely massive (> 100 Mg) metal- free stars formed 
fro m primordial gas ("population III" stars) may be an issue 
fe.g.lSchneider et alj|200a lYoshida et ai1l2004l ; iBaugh et ail 
120051 ). If so, using such an IMF in the modelling would lead 
to a reduction of the derived stellar masses and SFRs. 

We adopted iBruzual fc Charlotl (|2003h models with a 
metallicity of one-fifth solar and, for consistency, a Small 
Magellanic Cloud-type (SMC-type) d ust extinction law 
ijPrevot et al .111 9841 ; iBouchet et al]ll985h . Observational con- 
straints on the metallicity of z ~ 5 galaxies are scarce; how- 
ever, the metallicities derived from absorption lines in the 
spectra of eig th z ~ 5 LBGs in the Subaru Deep Field are 
Z ~ 0.2 Zq l)Ando et al.ll2004h . suggesting that our choice 
of 0.2 Zq is reasonably representative. For each individ- 
ual galaxy, best-fit parameters were derived from minimi- 
sation of the reduced \ 2 statistic. All model fitting followed 
standard procedu res applied in similar studies of high red- 
shift galaxies fe-g-IPapovich et al.ll200ll,l2006l:IShaplev etaD 
I200ll . l2005l : IForster Schreiber et al.ll2004l : lYan et alJl2006al ). 
We assumed a constant marginalisation for each of the mod- 
els. While our results are marginalised over the extinction 
law, metallicity, star formation history and IMF, we have 
actually explored a range of all of these parameters (see Ap- 
pendix rX| but the results presented correspond to our pre- 
ferred single values of these properties as described above. 
Detailed aspects of the modelling procedure will be de- 
scribed in a forthcoming paper (N.M. Forster Schreiber et 
al., in preparation, Paper II). 

We note that it is difficult to differentiate between the 
exact metallicity, extinction by dust and star formation 
histories for each individual galaxy based solely upon the 
goodness-of-fit because of degeneracies among these model 
parameters. Therefore, to determine the confidence intervals 
for the modelled properties, we ran 500 Monte Carlo sim- 
ulations for each object. We applied the same best-fitting 
procedure after perturbing the input broadband photome- 
try assuming the photometric uncertainties are Gaussian as 
indicated by our background noise fluctuations analysis (see 
Section[2j. The results of these simulations provide the prob- 
ability distribution in parameter space for each source. By 
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Figure 3. Composite probability distributions of age and (a) 
stellar mass, (b) star formation rate and (c) extinction for the 
robust sample of z ~ 5 LBGs as determined for each galaxy from 
our 500 Monte Carlo simulations. The points overlaid show the 
best-fit properties for each object in the sample and the large 
cross indicates the formal best-fit median properties of the entire 
ensemble. The overlaid contours indicate the 68% (solid line) and 
95% (dashed line) integrated probabilities on the ensemble prop- 
erties measured from the centroid of the distribution. In other 
words, 68% (and 95%) of all sources should have ages and stellar 
masses that fall within these regions. The bulk of the composite 
probability distribution lies at ages less than 100 Myr, as do the 
best-fit ages for the majority (two-thirds) of the sources. 
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combining these individual probability distributions, we de- 
rived __th£se_fc£_th£prorjerties of the ensemble of sources (see, 
e.g.. IPapovich et aDl200ll . for an analogous approach). Ex- 
ample combined probability distributions for the age, stellar 
mass, star formation rate and extinction for all of the galax- 
ies in the robust sample are shown in Figure [3] These dia- 
grams show the probability that any galaxy satisfying the 
Lyman-break criteria has of having those properties. From 
these we can therefore determine the characteristic proper- 
ties of the ensemble of galaxies without relying solely upon 
the mean of the (degenerate) best-fit solutions. 



4 ENSEMBLE PROPERTIES OF THE ROBUST 
SAMPLE OF z~5 LBGS 

4.1 Redshifts: Photometric and spectroscopic 

The LBGs in the robust sample have best-fitting photomet- 
ric redshifts in the range 4.60 to 5.54 (<z> me£ K an = 4.8) 
as expected given our initial photometric selection^. Our 
photometric redshifts agree exceptionally well (within 2a, 
or better), for both the high redshift LBGs and low-redshift 
interlopers, with spectroscopic redshifts where available (see 
Figure Hand Tabled]). 

As part of the public spectroscopic 
surveys of the CDF S conducted by ESO 
llVanzella et all 120051 . l2006l . 



expected for young stellar populations with no more than 
moderate dust obscuration. The results of our stellar evo- 
lutionary synthesis modelling with the parameters specified 
in Section [3] yields the following typical properties for the 
galaxies. These are typically strongly star forming galaxies 
with a median star formation rate of 40 Mq yr _1 . The me- 
dian photometric stellar mass is 2xl0 9 Mq, about a factor 
ten lower than Lyman-brea k galaxies at z ~ 3 |Shaplev et al.l 
[ioOll : IPapovich et al .1 f200 lh . The stellar mass estimates are 
the most robust of all derived properties, varying by a factor 
of 2-5 depending on the model assumptions (see Appendix 
lA"]l . While a large spread in the best-fit age is seen amongst 
individual sources in the sample the median age is relatively 
young, ~25Myr. The best-fit models imply only moderate 
extinction with a median value in the V-band of Av ~ 0.3 
mag. Models including extinction at this level fit the SEDs 
of the z ~ 5 LBGs better than models without extinction. 

As discussed in Section [3] we use the results of our 
Monte-Carlo simulations to characterise the properties of 
the ensemble of z~5 LBGs. Figure [3] shows the combined 
probability distributions for the parameters of mass, star 
formation rate and extinction as a function of age of the 
stellar population. While the ensemble probability distribu- 
tion in age extends over 1 Myr to 1 Gyr, there is a clear 
concentration at young ages. As indicated by the contours 
shown in Figure more than 68 per cent of the age dis- 
et all 120051 . I2006I . but also see vl.2 release tribution lies at ages less than 100 Myr, consistent with the 
|http: / /www.eso.org/science /goods /spectroscopy/CDFS-Maste rcaff^jority (two-thirds) of the galaxies having formal best-fit 
spectroscopy has been performed for 7 LBGs from our ages younger than this. In addition, Figure a) shows that 



robust sample, 6 of which are confirmed to be at z ~ 5. The 
low signal-to-noise ratio of the spectrum of the remaining 
robust LBG prevented assignment of a redshift. This high 
confirmation rate reinforces that our selection of z ~ 5 LBGs 
is robust. Spectroscopic redshifts were also secured for 17 
additional sources: 4 are stars, 4 are low-redshift interlopers 
(all agreeing with our photometric/morphological classifi- 
cation) and 9 are confirmed to lie at z ~ 5 and are part of 
our 'insufficient photometry' sub-sample which we expect 
includes genuine z ~ 5 LBGs. The right panel of Figure 
[J] shows the z-band magnitude and redshift distributions 
of the spectroscopically confirmed redshift 5 LBGs from 
the robust sample and those flagged as having insufficient 
photometry. The latter are slightly fainter than the robust 
sample while the redshift distribution is similar indicating 
bona fide fainter LBGs are among the sources flagged as 
having insufficient photometry. 



4.2 Stellar Masses, ages, star formation rates and 
extinction 

The rest-frame SEDs of the majority of our prime candi- 
dates are very blue from the far-UV to the visible, as is 

2 Using standard procedures l lMadau et al Jl996h we estimate the 
redshift range our selection criteria are sensiti ve to. This range is 
determ ined following the implementation in iLehnert fc Bremej 
(2003). Standard local galaxy templates were modified to ac- 
count for IGM opac ity along the line of sight short-ward of Lya 
ijMadau et al]|l996h . They were then further modified using the 
re dshift dependent Gunn-Peterson absorption of Lya as described 
111 iFan et al.ll|2005h . This results in star-forming galaxies lying be- 
tween z~4.6 and z~6 being able to satisfy our selection criteria. 



these young sources have stellar masses <10 Mq and typ- 
ically ~10 9 Mq, and Figure E|b) suggests they are the most 
strongly star- forming systems within our sample. This is un- 
surprising as they must have the lowest UV mass-to-light 
ratios in our sample given their low mass. 

The remaining third of our sample have best-fit masses 
of order 1O 1O M0 and ages older than a few hundred million 



(EEami et al. 2005; Schacrcr & Pello 2005; 


Mobasher et al.l 120051: 


Charv et al.l 120051: 


1: iDow-Hveelund et al. 


20051: Yan et al. 



7) than ours, where the presence of a discontinuity between 
the near and mid-infrared bands in the SED is identified 
with the rest-frame Balmer/4000A break (Figure [TJ. This 
break is indicative of more evolved stellar populations, where 
the emission from stars with ages of several hundred million 
years dominates over that from the short-lived massive stars 
that produce the UV emission. The absence or weakness of a 
discontinuity between the JG-band and 3.6/xm photometry 
points in the SEDs of the majority of our sample of galax- 
ies with younger best-fit ages strongly limits the possible 
contribution to the integrated stellar light by such similarly 
evolved stellar populations. Indeed, the lack of this disconti- 
nuity constrains their best-fit ages to less than ~100 million 
years. While old and massive systems are present in the 
sample, our results clearly indicate that a substantial frac- 
tion (>two-thirds) of galaxies at z~5 satisfying the Lyman 
break selection criteria are dominated by a young, intensely 
star-forming component. Similarly young ages and moder- 
ate masses (ages <45Myr and stellar masses ~1O 9 M0) are 
found for IRAC undetected z ~ 6 ^- band-dropouts (z ~ 6) 
LBG candidates recently reported by I Yan et ahl (|2006al ). 
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Figure 4. Left: This figure shows the agreement between spectroscopic and photometric redshifts of all objects which satisfy our 
selection criteria which are confirmed to be galaxies. All filled symbols are for z~5 LBGs from our robust sample. Two a error bars on 
the photometric redshift are overplotted. The figure excludes a low redshift interloper which has an x-ray counterpart (XCDFS265) and 
is likely to be AGN-hosting, therefore the results of the stellar evolutionary synthesis modelling does not make sense in this case. As a 
guide, the line z sp =z p f l is overplotted. Right: z-band magnitude distribution of the spectroscopically confirmed LBGs. As in Figure [2] 
the black/shaded histogram corresponds to galaxies in the 'robust' sample. 



Table 1. Breakdown of photometric and spectroscopic redshifts for all of the LBG candidates. The z v hot column refers to the nominal 
best-fit redshift determined from our modelling and the z spec to spectroscopic redshifts from the literature. 







z phot 






Zspec 








Total 


z<4 


z>4 


Attempted 


z<4 


z>4 


Unassigned 




(la) 


(lb) 


(lc) 


(2a) 


(2b) 


(2c) 


(2d) 


UV Selection 


109 


23 


86 


34 


8 


15 


11 


Robust sample 


21 





21 


7 





6 


1 


Insufficient Photometry 


17 


1 


16 


11 





9 


5 


Star/QSO 


19 





19+ 


5 


4 





1 


Interlopers 


22 


22 





8 


4 





4 



Columns l(a-c) refer to the results of the SED modelling for all candidates (109) that satisfy our original UV-selection, and columns 
2(a-d) for those of the 109 candidates which have been spectroscopically follow-up. 
(la) Total number of sources satisfying our initial selection criterion. 

(lb) & (lc) The numbers of LBG candidates with photometric redshifts above and below z=4. 
(2a) Total number of LBG candidates from our initial sample that have had spectroscopic measurements. 
(2b) Number of candidates for which the spectra were of insufficient quality to definitively assign a redshift. 
(2c) & (2d) The numbers of sources with confirmed spectroscopic redshifts above and below z=4. 

+ Best-fit redshifts of z p ^ ot >4 are produced for all objects classed as stars/QSOs due to the the strong break in their SEDs being 
identified with the Lyman-break and because their SEDs are fit with inappropriate star-forming galaxy templates. 



In this section we have reported on the typical physical 
properties for the sample of luminous z ~ 5 LBGs for our 
adopted modelling assumptions. We have extensively inves- 
tigated the effects of varying the input model parameters on 
the derived properties and find the key properties of young 
ages and moderate masses are robust under a wide range 
of model assumptions. We discuss the effects of varying the 
input parameters in the Appendix (see also N.M. Forster 
Schreiber et al. in preparation). 



5 COMPARISON TO LBGS AT REDSHIFT 3 



5.1 Stellar Masses, Ages, Star-formation Rates 
and Extinction 

In this section we compare the properties of our z ~ 5 sam- 
ple with thos e pub lished for the z ~ 3 LBG sample from 
IShaplev et al.l (|200ll ) which is matched in rest-frame UV lu- 
minosity. The SEDs for both samples span the same rest- 
frame wavelength range and modelling shows significant dif- 
ferences in the stellar masses and ages derived even for 
identical modelling assumptions. Because of this and the 
matched selection, we can only ascribe these differences to 
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Figur e 5. Histograms (peak normalised) showing the distribution of best-fit properties of 74 R<25.5 z~3 LBGs from lShaplev et al,l 
(2001) (shaded grey) together with the best-fit properties of our robust sample of z~5 LBGs (line filled histograms). These properties 
correspond to iden tical model assumpt i ons: c onstant star formation, 0.1-100 Mq Salpeter IMF, Calzetti extinction law and solar mctal- 
licity templates of iBruzual fc C hariot ( 20031) . The vertical lines indicate the median values for the z ~ 3 sample (dashed line) and the 
z~5 sample (solid line). LBGs at z~5 are typically younger (<100Myr) and are less massive (fewxl0 9 MQ) than z~3 counterparts of 
similar luminosity (see text for a detailed discussion) . 



intrinsic differences in the properties of LBGs at these two 
epochs. 

The i-band selection limit of us < 26.3 chosen for our 
z ~ 5 sample is matched to the typical magnitude limit of 
spectroscopically confirmed samples of LBGs at z ~ 3. This 
limit (TZab <25.5), is sensitive to z ~ 3 LBGs that are ~ 1 
mag fainter than m* z ~ 3 = =24.48 (Stcidel c t~ai1 1 19991 ). 

1700A.AS ' ' 

The 1Z— band maps to 1700A in the rest-frame at z ~ 3. This 
enables us to compare the ensemble properties of our robust 
sample of z ~ 5 LBGs to those obtained for 74 z ~ 3 LBGs 
with apparent m agnitudes grea t er tha n 0.1m* from the sam- 
ple analysed by IShaplev et all l|200ll ). As well as matching 
our samples by magnitude, in order to make a fair compar- 
ison we have modelled the SEDs of our ro bust sample of 
z ~ 5 SEDs using the same assumptions as IShaplev et all 
Specifically, we used solar metall i city m odel spectra gener- 
ated with the lBruzual fc CharlotJ (|2003l ) population synthe- 
sis code, a constant star formation history, a Salpeter IMF 



and a Calzetti attenuation law ()Calzetti et al.ll2000l ) . Note 
that this is different a extinction law and metallicity than 
those used in previous sections, but allows a direct compar- 
ison between the two samples. Any differences in the results 
indicate intrinsic differences in the properties of the two sam- 
ples, either in properties like mass and star formation rate, 
or that different extinction laws or metallicities are required 
at the two epochs. 

Figure [3] shows the result of this comparison for the 
best-fit stellar masses, ages, extinctions and star formation 
rates obtained for a con stant star formation history from 
the z~3 LBG sample of IShaplev et al.l (|200ll ) (shaded his- 
tograms) and the results for our robust z ~ 5 sample (line 
filled histograms). While the derived extinctions are well 

3 They assume a slightly higher upper mass cut-off of 125 Mq 
than our assumed 100 Mq , however this slight increase has an 
insignificant effect on the generated models and resultant prop- 
erties. 
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log I0 ( SFR [M yr" 1 ] ) A v [mag] 

Figure 6. The shading and lines are as in Figure \E\ but now showing the distributions of properties derived from our Monte-Carlo 
simulations of the properties of z ~ 5 LBGs from our robust sample (line filled histograms) obtained with 0.2 Zq templates and a SMC 
extinction law (see text for details). 




log 10 ( Stellar Mass [M ] ) log 10 ( SFR [M yr _1 ] ) 

Figure 7. This figure shows that the distribution of properties derived from our Monte Carlo simulations (line filled histograms) trace 
the distributions of nominal best-fit (shaded histogram) properties well. Stellar masses and star formation rates are shown, ages and 
extinctions show similar agreement. The solid vertical line indicates the median value. 
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matched, this comparison clearly shows significant differ- 
ences in the ages, masses and star formation rates between 
the samples at z ~ 3 and z ~ 5. The z ~ 5 sample shows 
young typical ages (<10Myr), in contrast the z~3 LBGs 
are significantly older (~320Myr). The typical stellar mass 
of z~5 LBGs is ~5-10 times lower than the z~3 sample 
(few xlO 9 Mq). Extremely high star formation rates are es- 
timated (~500 M© yr" 1 ) for the z ~ 5 LBGs, a factor of ten 
higher than the z ~ 3 sample. Since we have compared the 
properties of similarly selected galaxies, based on the same 
rest-frame UV measurements, to the same UV luminosity, 
which have been modelled under the same assumptions, 
these differences can unequivocally be ascribed to inherent 
differences in the properties of the two samples of LBGs. 

Moreover, this comparison demonstrates that adopting 
Calzetti extinction law and solar metallicity may not be ap- 
propriate for this sample of bright z ~ 5 LBGs. Implausibly 
young ages (<10Myr) and excessively high star formation 
rates (>lO 3 M0yr _1 ) are seen for a large fraction of the 
z ~ 5 population. The fact that the ages are so signficantly 
younger than the z ~ 3 LBG population (see top left panel 
Figure [5| is an indication that the SEDs of z ~ 5 LBGs are 
even 'bluer' than their z ~ 3 counterparts. The combination 
of the adopted extinction law and the high metallicity tem- 
plates results in these extremely young ages and high star 
formation rates (see the discussion in Appendix 1X1 and also 
IPaoovich et~aDl200ll L 

As discussed in Section [3] evidence suggests our consis- 
tent choice of 0.2 Zq and the SMC extinction law is suitable 
to mo del the blue SEDs of our z ~ 5 LBGs. IShaplev et alj 
l|200ll ) have discussed the suitability of the Calzetti plus so- 
lar metallicity models for z~3 LBGs. In Figure [6] we show 
the comparison between the properties of LBGs at z ~ 3 and 
z ~ 5 under these appropriate (but different) assumptions. 
We have elected to show the probability distribution of prop- 
erties obtained from our Monte Carlo simulations of the z ~ 5 
LBGs to give less importance to the formal best-fit solution 
which can be sensitive to the modelling assumptions. Since 
the probability histograms follow the same distribution as 
the best-fit parameters themselves (see example for stellar 
mass and SFRs in Figure 0, our comparison using the prob- 
ability distributions is reasonable. 

Figure [H] shows that the distributions of star formation 
rates are extremely well matched, however, there are clear 
differences in the distributions of stellar mass, age and ex- 
tinction. The typical stellar mass of LBGs is a factor of ten 
lower at redshift 5 than at redshift 3. Similarly, the dis- 
tribution of ages shows more sources have ages <100Myr 
in the z ~ 5 sample than the z ~ 3 sample. There is an 
apparent difference in the distribution of extinctions but 
this can largely be attributed to differences in the mod- 
elling assumptions (Appendix \K§ . This is corroborated by 
the fact that the extinction distributions overlap almost en- 
tirely when the same modelling assumptions have been made 
as we have shown in Figure [3] Furthermore, a comparison 
to the results for sim ilarly luminous z ~ 3 LBGs studied in 
IPapovich et alj (|200ll ffl assuming a Calzetti extinction law 



4 iPapo vich ct a lj model z ~ 3 LBGs with the same metallicity 
l|Bruzual fc Charlod Ejoolh templates and IMF as we have as- 
sumed. However, they give best fit results for individual galax- 



but with the 0.2 Zq templates, shows similar extinction for 
the z ~ 3 LBGs as is seen in Figure [H] Thus the difference in 
the extinction between z ~ 5 and z ~ 3 LBGs seen in Figure 
[B]is mostly a consequence of the adopted extinction law and 
not the metallicity of the templates. 

IPapovich et alj (|200ll ) have comprehensively discussed 
the effects of different model assumptions on the derived 
properties of a faint sample of z ~ 3 LBGs. Their combined 
probability distributions highlight the problems associated 
with using the Calzetti extinct ion law with the low metal- 
licity iBruzual fc Charlotl < |2003h templates. In their Figure 
9b, they show the probability distribution resulting from 
the SED modelling assuming a Calzetti extinction law with 
the bluer 0.2 Zq metallicity templates. A clear extension of 
sources with extremely young ages (< lOMyr) and high 
Av (~2) is seen. This extension almost disappears under 
a Calzetti law with Zq, and with the SMC law for either 
the 0.2 Zq or 1 Zq models (compare their figure 9b with 9d, 
12b and 12d for a Salpeter IMF). 

Thus, adopting the Calzetti law with the solar or low 
metallicity templates would only act to reduce the derived 
ages of the z ~ 5 LBGs and results in stellar masses that are 
5-10 times lower than z ~ 3 LBGs, i.e. leaving our primary 
conclusions of young ages and moderate masses for z ~ 5 
LBGs unaffected. Moreover, under all suitable assumptions 
for the extinction law and consistent metallicities (i.e. 1 Zq 
with the Calzetti law, 0.2 or IZq with SMC), we find that 
galaxies with ages (<100Myr) and masses (~ 10 9 Mq) dom- 
inate samples of bright LBGs at z ~ 5. Our analysis of the 
probability histograms and distributions shown in Figures 
[3] and [6] shows that this fraction to be >70 per cent. How- 
ever, while z ~ 3 LBGs with such yo ung ages and similar 
stellar masses are n ot unknown (e.g. IPapovich et alj|200ll ; 
Shapl ev et alJl200ll ). they are less than half as prevalent in 
LBG samples at z~3 than at z~5 (<30 per cent). This 
unequivocal change in the age and mass of the dominant 
population in luminous LBG samples separated by ~1 Gyr, 
implies that at z ~ 5 is an epoch where the majority of LBGs 
are recently formed and are in the process of accumulating 
their first significant mass. 



5.2 Duty cycle and bias 

We have shown in the previous sections that the dominant 
population of LBGs in our z ~ 5 sample is younger and of 
lower mass than that in z ~ 3 LBG samples with similar 
UV luminosities. Since the majority of LBGs at z ~ 5 are so 
young, it is inevitable their detectability is highly stochastic 
and this stochasticity impacts upon the interpretation of 
other statistical properties such as their clustering strength, 
number density and bias. 

ies only with the Calzetti extinction law. Since Papovich et al. 
give best fit results for a range of star formation histories we con- 
sider only those fit with constant star formation and those with 
best-fit ages less than the e-folding time-scales of their 'best-fit' 
exponentially decaying star formation history model. The latter 
condition for age<T gives similar results to constant star forma- 
tion. In addition, we considered only those LBGs from this faint 
sample which are matc hed in luminosit y to the z~5 samples, 
V<24.8. Ten LBGs from [Papovich et alj satisfied these criteria. 
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iLehnert fc Bremerl (|2003T ) showed that the co-moving 
number density of z ~ 5 LBGs selected in a similar manner 
to our robust sample was a factor of three lo wer than that 
of sim ilar luminosity LBGs at z ~ 3 and z ~ 4 |Steidel et al.l 
1999). If the LBGs at both epochs had identical properties 
then this can be directly related to evolution in the comov- 
ing number densities of the halos that they occupy. Clearly, 
this naive picture is complicated by the clear difference in 
properties of the two samples in the previous section. More- 
over, the younger ages of the higher redshift sample imply 
a shorter duty cycle for the LBG, related to the length of 
time a galaxy would be detected as a UV luminous LBG. 

Given that the average lifetime of the z ~ 5 LBGs of 
^25Myr and that the lookback time spanned by our sur- 
vey (325 Myr over the redshift range of the z ~ 5 LBGs 
4.4 < z < 5.6), suggests that luminous members of the z ~ 5 
LBG populations have a duty cycle of . 08. A similar es- 
timate for z~3 LBGs by IShaplev et all (1200 il l, renders a 
duty cycle of 0.5, significantly longer than at z ~ 5. So across 
the two sample volumes, the higher redshift sample rep- 
resents systems ^thirteen times as numerous as those de- 
tected, whereas the lower redshift sample misses only half of 
the potential LBGs in the volume. Clearly, the young ages 
of bulk of the high redshift sample have a significant effect 
on the bias of these sources. 

LBGs have been long kn own to be a signifi- 



cantly clustered population (e.g. lAdelberger et all 1199 



IGiavalisco et al.ll 19981 : ["Ouchi et al.ll2004bl : ILee et al.ll2006l 

Clustering analysis of 23.5 < IZab < 25.5 LBGs at z ~ 3 
suggests that L BGs reside in dark ma tter halos of masses 
lO 11,2 " 11 ' 8 Mq (jAdelbereer et all lioolih consistent with re- 
cent estimates of the dynamical virial mass of the DM 
halos of LBGs determined rotation curves of z ~ 3 LBGs 
l| Weatherlev fc Warrenll2005l : iNesvadba et alllioOrih . For z~ 
5 LBGs. ILee et al.l l|2006h suggest slightly lower halo masses 
may be plausible {M^aio ~ 10 11 Mq) however the uncer- 
tainties are too large to see a significant difference between 
the z ~ 3 and z ~ 5 popul a tions. This conclusion is simi- 
lar to that of lOuchi et ail (|2004tJl who found an increase 
(of marginal significance) in the correlation length with in- 
creasing redshift which is consistent with a constant halo 
mass of ~10 12 M over z=3-to-5. Again, however, the un- 
certainties, especially for the highest redshift LBGs, were 
large. Given this, and in the absence of dynamical mea- 
surements of z ~ 5 LBGs, we adopt t he z ~ 3 mass range 
( 10 ii.2-ii. 8 M m . lAdelberger" et al . 2005) for the dark matter 
halos for z ~ 5 LBGs. Given that the comoving number den- 
sity of underlying dark matter halos of sufficien t mass to host 
such LBGs is evolving strongly from z = 5 to 3 l|Mo fc White! 
l2002h. the bias of th e LBGs must be increasing with redshift 
jOuchi et aill2004bl ). 

The uncertainties in the clustering statistics of LBGs 
at z ~ 5 are still relatively large. The determination of ac- 
curate halo masses at z ~ 5 will come from a analysis of 
the clustering statistics of a sufficiently large sample of ei- 
ther spectr oscopically-confirme d or photometrically-robust 
LBGs (e.g. lOuchi et al.ll2004bh . As we have shown in this 



5 Further analysis on the clustering strength of z ~ 5 LBGs will be 
presented in a forthcoming paper A. Verma et al. in preparation, 
Paper III 
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Figure 8. Evolution of the star formation rate density with red- 
shift. The SFRD determined from this study is shown (filled 
black circle, see text for details of its estimation) with respect to 
previous estimations and theoretical models. We emphasise the 
high redshift tail of the SFRD(z) and therefore only a handful 
of the numerous measurements at z<2 are shown. For consis- 
tency, all data and models are normalised to a Chabrier/Kroupa 
IMF. All data are corre spond to an i ntegra tion of the z~3 
LB G luminosity function Steid el et al ] dl999h to 0.1 L*, except 
for I Giavaliseo et al. | a 2004aF and lBouwenset al. I ll2005li which are 
to 0.2L* and 0.04L Z — 3*, respectively. Measurements from the 
literature were estimated using scaling relations to the UV lu- 
minosity and account for an extinction of E(B-V)=0.15, as per 
iNaeamine et al.l a 2006t) . 



paper, the latter requires a sample defined using both opti- 
cal and infrared bands to exclude low redshift galaxies. Until 
this is done comprehensively, clustering statistics from sam- 
ples selected by optical photometry alone or those that only 
include shallow infrared data, should be interpreted with 
caution. 



5.3 Evolution of the global star formation rate 
and stellar mass densities 

5. 3. 1 Global star formation rate density 

The evolution in the star formation rate density (SFRD) 
with redshift has been a crucial tool in our understand- 
ing of the evolution of galaxies in an ageing Univ erse. The 
ground-breaking 'Madau-plot' ijMadau et al. Il996h and ver- 
sions thereafter have clearly defined the era of peak star for- 
mation activity in the Un iverse (z~ 1-2, e.g. I Connolly et ah! 
1 19971 ; iHeavens et al.ll2004h . however constraining the SFRD 
at higher redshifts is still an issue of debate with direct 
observed estimates suggesting a decline for z > 3, while 
estimates based on integrating ov er the luminosity func- 
tion and theoretical predictions ijSomerville et al . 2004; 
iNagamine et al.l I2006T ) imply a consta nt value or shallow 
decline beyond the peak at z ~ 2 (e.g. Steidel et all 19991; 
Bunker et all 120041: lOuchi et all l2004al: IGiavalisco et all 



2004al ; iBouwens et al 



2005 



also iHopkins fc Beacoml 



Bouwens fc Illingworthl 2006: 



20061 f or a compilation and 



refere nces therein). The model of iHernauist fc Springell 
(2003) and the TVD hydrodynamical simulations of 
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Naga mine et alJ i|2006l ) suggest the SFRD might still be in- 
creasing at z ~ 5. 

Our results from the SED modelling of z ~ 5 LBGs has 
an impact on this. The SFRD derived from the best-fit 
star formation rates for LBGs in our sample is compared 
in Figure [8] to several recent observational determinations 
and model predictions from the literature. This is deter- 
mined from the total SFR for our robust sample, corrected 
for the number of true z ~ 5 LBGs we expect in our to- 
tal sample of 109 sources (all except the 40 per cent in- 
terloper fraction), incompleteness at our magnitude limit 
(90 per cenlQ) and accounting for galaxies to 0.1 L^Q over 
our survey volume. Nominally, the survey volume is given 
by the area of the field, multiplied by the depth in red- 
shift space between the minimum and maximum redshifts 
that can satisfy our initial selection criteria. However, not 
all sources would remain in our magnitude limited sample 
because of cosmological dimming and increased intervenin g 
hydrogen opacity with redshift l|Lehnert fc Bremerl 120031) . 
Making the assumption that all sources are at the maxi- 
mum possible redshift to remain in the sample (i.e. they are 
close to the magnitude limit of iab = 26.3) gives a volume 
probed of ~40 per cent the maximum notional volume (con- 
sistent with fractions det ermined from detailed calculations 
iLehnert fc Bremerl [2003. and E. Stanway, private commu- 
nication). The star formation rate density accounting for 
this effective survey volume is shown if Figure [H] For consis- 
tency with the models, the SFRs determined have been cor- 
rected to a Chabrier/Kroupa IMF, rather than the Salpeter 
IMF we ha ve used in the modell ing, by dividing by a fac- 
tor of 1.6 (|Nagamine et alj|2006h . The high SFRD deter- 
mined supports no more than shallow decline in the SFRD 
at high redshift, consisten t wit h the theoretical predictions 
of INagamine etafl (|2006l ) and iHernquist fc Springell 120031 . 
however lies above the semi-ana lytic accelerated quiescent 
model of ISomerville et al.1 (120011 ). 

Star formation rates can also be derived by scaling the 
intrinsic rest-frame UV luminosities of sources. The addi- 
tional data points shown in Figure [5] are derived as such. 
For starbursts younger than ~ 100 Myr, the ratio of star 
formation rate to UV luminosity declines with age, level- 
ling off at ~ 100 Myr. This is due to the time it takes to 
build up a UV luminous stellar population in such a galaxy, 
the levelling-off occurring once their birth and death rates 
match. In part owing to this effect, az~5 sample matched in 
intrinsic UV luminosity with az~3 sample will have higher 
star formation rates for the typically younger ages we infer. 
Here we note that the UV luminosity d ensity declines by a 
facto r of ~ 3 between z = 4 and z = 5 ([Lehnert fc Bremerl 
2003). Because of the younger ages at higher redshifts, this 
must be an upper limit to the decline in SFRD, at least for 
the more luminous LBGs. A detailed exploration of this is- 
sue is left to a forthcoming paper (N.M. Forster Schreiber 
et al., in preparation, Paper II). 



6 Gwyn, determined from the vl.O release of the ACS data see 
http: / /www, astro, uvicca/gr ads /gwyn/virm os/cdfs/index 5 , html 
' assuming the z ~ 3 luminosity function of Stcidel et al. (1999) 
(a=-1.6) since the faint end slope of the z ~ 5 luminosity function 
remains uncertain. 
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Figure 9. Evolution of the stellar mass density with redshift. The 
left y-axis is plotted in the dimensionless Q*=p* / p c * parameter 
and the right shows the percentage of stellar mass that has been 
assembled i n the Universe with respect to the z=0 estimate from 
ICole et al.l j200ll) (open star). All data and models correspond to 
a Chabrier/Kroupa IMF. The upper and lower limits determined 
from this study are plotted as (filled circles). The models are the 
integral of the SFRD(z) which are modified to account for mass 
loss from evolved stars. 



5.S.2 Global stellar mass density 

The evolution in the global stellar mass density with red- 
shift places strong constra ints on models of the st ellar mass 
assembly in the Universe ([Dickinson et al ] |2003bh . Approx- 
imately 50 per cent of the stellar mass in galaxies had 
formed by z^ 1-1.5, with ~10 per cent in place by z ~ 3 
[jPickinson et al.1 l2003bl : iRudnick et all 12001 120061 ). Accu- 
rately constraining the evolution of this property with red- 
shift requires mass-selected (or limited) samples at differ- 
ent epochs. While our sample of luminous LBGs is not 
mass selected, we can nevertheless place limits on the for- 
mation of stellar mass soon after the end of reionisation. 
We have determined the stellar mass in our survey vol- 
ume after applying the incompleteness corrections described 
in the preceding section (with the exception of the correc- 
tion to 0.1L*). The stellar mass density at z~5 is ~ 
1.5 x 10 6 M Q Mpc~ 3 (Q* ~ 1.1 x 10~ 5 ), for a Salpeter 
IMF. For consistency with the published modelling, e.g. 
INagamine et al.l l|2006h . we reduce this value by 1.4 to ac- 
count for the difference between the Salpeter IMF and the 
Chabrier/Kroupa IMF used in the models, as shown in Fig- 
ure[9] This limit represents ~0.3 per cent of the total stellar 
mass of galaxies seen today and it is consisten t with recent 
estima t ions at high redsh ift Istark et al.l (|2006l ); lEvles et al.l 
l|2006l ): lYan et al.1 l|2006al ). 

The stellar mass density derived is a lower limit partly 
because of the stochasticity discussed in Section 15.21 The 
relatively short duty cycle of these sources suggests that up 
to ~12-in-13 LBGs at z~5 may be missed in our survey. 
Accounting for this incompleteness places an upper limit on 
stellar mass density at this epoch for our sources (see Figure 
[9|. Even after accounting for this, our derived upper limit 
falls below the predictions of the evolution of mass density 
with redshift obtained from the integration of the models 
that plausibly explain the SFRD(z) to high redshift (see 
Figures |5J). This is to be expected given our sample is not 
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selected to be mass complete. Extending the mass function 
to lower masses would again increase the mass density but at 
a lower level (a factor of a few) compared to the duty cycle 
correction. Moreover, at z~2 as much of the st ellar mass den- 
sity contained within UV-selected samples l|Shaplev et al.l 
120051) lies in the pop ulation of massive distant red galaxies 
l|Rudnick et~aLll2006h , which are largely missed by UV selec- 
tion techniques. Comparably massive galaxi e s (>1 11 M©) 
in the UV selected sample of IShaplev et al.l (|2005T j are an 
order of magnitude less numerous than the DRGs. Whether 
such massive galaxies dominate the mass density at z ~ 5 
remains to be established. 

The apparent discrepancy between these theoretical 
predictions and the observed estimates for z > 1 has been 
attributed to incompl eteness of the observational data (e.g. 
iNagamine et al.l l2006h . or deficiencies in the models, pre- 
dicting too much stellar mass at high r edshif t, despite 
reproducing the SFRD(z). iRudnick et al.l (|2006l 'l compre- 
hensively discuss both observational estimates and model 
predictions highlighting the potential causes for the dis- 
crep ancies from both perspe ctives (see also the discussion 
in iGwvn fc Hartwickl |20 05). Notwithstanding, based on 
the observational data alone, our upper limit implies that 
~1.5 per cent of the present-day stellar mass density could 
have already been been assembled in luminous LBGs in 
the fi rst 1.2 Gyr, rising b y an order of magnitude ~1 Gyr 
later l|Rudnick et al., ' 2006). This suggests that once the first 
galaxies had formed, star formation and stellar mass assem- 
bly proceeded at a prodigious rate in the early Universe. 

5.4 Mass and luminosity densities 

We measured the sizes of LBGs in the robust sample us- 
ing the ACS z-band images. Using the spectroscopic and 
photometric redshifts, the mean (proper) size of the z ~ 5 
LBGs is found to be 1 kpc. This is a 40 per cent decrease 
from the mean size of LBGs at redshift 3 (< r e >~ 1.7 kpc), 
which is fully consist ent with the known s i ze evolution of 
UV-selected galaxies ((Ferguson et al 1 l2004l : iBouwens et al.l 
120041: iBremer et al.ll2004h . 

iTruiillo et alj (|2004T ) found that there was no or only 
weak evolution of the stellar mass densities (mass/r^) of 
galaxies with masses > 10 10 Mr over the redshift range of 
to ~3. While the results of ITruiillo et al.l (12004 ) are in- 
complete for lower mass galaxies, we that the stellar mass 
de nsities for our LBGs at z=5 are comparable to that found 
bv lTruiillo et al.1 l|2004l ). However, their UV luminosity den- 
sities (luminosity /rj!) are significantly higher than LBGs at 
z ~ 3 with z ~ 5 LBGs producing as much luminosity as z ~ 3 
LBGs in one-fifth of the volume. This difference cannot be 
attributed to a low surface brightness, extended component 
since the ACS images are significa ntly deeper than our m ag- 
nitude limit for selection (see also lBouwens et al.ll2004h . 

The fact that the stellar mass densities are compara- 
ble to star-forming galaxies at low redshifts, while the ob- 
served luminosity densities are significantly higher, is indica- 
tive that the mas s-to-light ratio of su ch galaxies is evolv- 
ing with redshift l|Truiillo et al.ll2004h . LBGs (and DRGs) 
at z~2-to-3 have <M/Ly > ratio of ~0.5- 1 (|Truiillo et al.1 
l2004l : IShaplev et al.ll2005h while we find by z ~ 5 LBGs have 
<M/Lv> ratio of 0.1. In the hierarchical framework, stellar 
mass assembly is likely to have occurred in episodic bursts. 



Low mass-to-light ratios, such as those we are finding for 
the z ~ 5 LBGs, are expected for galaxies at high redshift 
which are forming stars in a recent burst. These findings are 
consistent with our contention that these are young (and 
likely without a significant intermediate age population, see 
Appendix [Bj, unextincted intensely star-forming galaxies. 



5.5 Lyman alpha emitting galaxies 

The spectra of 9 of 15 LBGs with spectroscopic redshifts 
show an emission feature with an asymmetric line profile 
and continuum break consistent with the being Lya. For 
z ~ 5 - 6 LBGs, lAndo et al.l l|2006h find that luminous LBGs 
are Lya deficient, with the separation between those which 
are Lya emitters occurring at >M* (M « ~-21.5 to 

J ° ~ v AB.1A0QA 

-21.0 mag). They postulate that more luminous LBGs are 
more chemically evolved systems where dust absorbs the 
Lya emission, although alternative explanations cannot be 
rul ed out. This is consisten t with the ideas put forward 
by iMori fc Umemural l|2006h who suggest that Lya emit- 
ting and Lyman-break galaxies are evolutionary stages in 
the chemical enrichment of galaxies that trace the evolu- 
tion of primeval irregulars to present day ellipticals. From 
N-body hydrodynamical simulations with stellar population 
synthesis they state that the Lya is produced by cooling 
shocks in the early, sub-solar metallicity evolutionary phase 
for ages < 300 Myr. Then, up to 1 Gyr, the stellar continuum 
increases and the metallicity becomes solar. 

We would therefore expect to see a trend in our de- 
rived ages and the equivalent width of Lya in our sam- 
ple with older galaxies having lower equivalent widths. Us- 
ing the ESO-FORS2 spectroscopy, we have crudely esti- 
mated the Lya equivalent widths (ignoring underestimation 
caused by the blue-shifted absorption due to intervening 
opacity along the line of sight) of the Lya emission to be 
-20 EW(Lya) ^ A for all LBGs in the robust sample, ir- 
respective of absolute magnitude (-24.0<M ■ <-20.5). 

^ a v AB, 1400A ' 

However, as spectroscopy has only been obtained for the 
brightest LBGs, it is not surprising that we do not see the 
extremely faint Lya emitting galaxies with high equivalent 
widths (<-40A) seen for M « >-21.5. There is no sig- 

nificant correlation between the strength of Lya equivalent 
width with the nominal best-fit ages. We do not find any 
strong difference in the distributions of mass or star forma- 
tion rate suggesting these properties do not influence the 
presence or absence of Lya. Furthermore, we find no evi- 
dence that the 'break-only' galaxies have a higher average 
extinction than the galaxies exhibiting Lya in their spectra 
and therefore attenuation by dust may not be the sole rea- 
son why there is a lack of Lya emission in the break-only 
LBGs. Rather local conditions must influence the presence 
and strength of Lya emission. It may well be that young 
galaxies have a high covering fraction of neutral gas which 
becomes successively more ionised as the galaxy ages. The 
effect of winds and the geometry of the gas may also play 
a crucial role in the range of ages seen in galaxies without 
significant Lya emission. 
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6 DISCUSSION 

We present detailed analysis of the properties of LBGs at 
z ~ 5 with robust and reliable photometry spanning the rest- 
frame UV to optical. Through a probability analysis we are 
able to constrain the properties of the ensemble and find 
that 'y°mig' (<100Myr) and moderately massive (~ 2 x 
10 9 M©) galaxies are in fact more than twice as prevalent at 
z ~ 5 (70 per cent) than at z ~ 3 (30 per cent). This increased 
fraction of young and moderately massive galaxies suggests 
that at redshift 5 we are seeing an era of widespread stellar 
mass assembly in the early evolution of galaxies. Based on 
these properties, we discuss three key characteristics of this 
population that emerge from our analysis. 

6.1 Progenitors of present day early type galaxies 
or bulges 

The luminous z ~ 5 LBGs studied here are undergoing a pe- 
riod of intense, recent star formation that dominates their 
rest-frame UV-to-visible emission. Their compact sizes, as 
measured from the HST images, show that the UV emit- 
ting regions are small (mean and median half-light radius 
of ~lkpc). In their short lifetimes, each of these galax- 
ies has already assembled >2x 10 9 Mq of stars within this 
radius. The inferred stellar mass surface density (/x* ~ 
6 x 10 8 Mq kpc -2 ) is comparable to the bulge and spheroidal 
components seen in prese nt-day typical (L* ) galaxies (e.g. 
< /i* >~5xl0 8 kpc~ 2 . [Gaiaz et alJliool ). 

From a study of nearly 400,000 low redshift galax- 
ies drawn from the Sloan Digital Sky Survey (SDSS), 
iKauffmann et all have shown that significant differences be- 
tween the structural parameters of early and late type galax- 
ies are highly dependent on their mass surface densities, 
rather than the total mass of the galaxies. Specifically, they 
find that the concentration parameter (C > 2.5) and stellar 
mass surface densities (/i* > 3 x 10 s Mq kpc -2 ) correspond 
to the regime of galaxy spheroids and bulges, independent 
of the total stellar mass of these systems. Above this thresh- 
old of mass surface density, they find that star formation 
is increasingly suppressed in local galaxies (the specific star 
formation is l ow) and must hav e cease d many Gyrs ago. As a 
consequence, IKauffmann et al. I (|2006l ) hypothesize that with 
increasing compactness and surface density of the galaxy, 
stars were formed in short, vigorous epis odes at high red- 
shift, with extended periods of inactivity (|Kauffmann et all 
120061 . characterized this as the gas consumption time scale, 
tcons oc fJ-Z 1 )- This analysis is particularly appropriate here 
since we do not find high stellar masses which is often taken 
as evidence for being progenitors of early type galaxies; we 
instead find very high mass surface densities. Given the cir- 
cumstantial evidence that the z ~ 5 LBGs in our sample 
drive winds (see Section [6.2. have young ages and short 
duty cycles, it is also plausible that these galaxies will grow 
in strong bursts of star-formation. So while we cannot know 
the evolutionary path of this population of galaxies, their 
overall mass surface densities suggest that they are the pro- 
genitors of early type galaxies or bulges, and favours the 
"inside-out" galaxy formation scenario. 

The star formation that dominates the UV-to-visible 
SEDs of the z ~ 5 LBGs in our sample is likely to be such 
a burst of star formation in which a significant fraction of 



the stellar mass of the system has been formed. A simple 
estimate of the relating the dynamical time-scale to the mass 
and size of a galaxy is 

_ r „ y^r^ 

l-dynamical — y ~ Ql/2 ' ' 

where V is the characteristic (virial) velocity of the sys- 
tem, rj is a geometrical correction factor that depends on the 
mass distribution within the system, r is the radius or physi- 
cal scale of the system, M is the mass of the system, and G is 
the gravitational constant. Depending upon the mass distri- 
bution for a given ga laxy, r\ is as 3-5, (Bi nnev fc Tremaind 
1 19871 : iRix et al] 119971 ) . With a median measured half-light 
radius of ~lkpc and a mass of 2 x 1O 9 M0, we estimate 
that the dynamical time-scale of the galaxies in our sam- 
ple is of order 20 Myr. Therefore, the intense star formation 
in the majority of the z~5 LBGs (age~100 Myr) has been 
typically proceeding for approximately one to ~5 dynam- 
ical time-scales. Starbursts, the youngest, large-scale star 
star forming e vents seen in nearby galaxies have compara- 
ble durations ( Lctmcrt & He ckmanl 1 19961 : iKennicuttl Il998l : 
iForster Schreiber et alj l2003). A system that has assembled 
a significant fraction of its stellar mass over the last few 
dynamical time-scales is in essence a galaxy in the process 
of formation. Thus, the young galaxies in our sample are 
likely to be members of the long sought-after population of 
primeval galaxies in the process of significant stellar mass 
build up. 

Further support of the progenitor scenario arises from 
analysis of the clustering of z ~ 5 LBGs. The clustering 
length of z ~ 5 LBGs implies that L>L* LBGs at z ~ 5 reside 
in dark matter halos with masses of order 10 11 " 12 M 
ijOuchi et al.ll2004bl : iLee et al.ll2006h . While this is compa- 
rable to the dark matter mas s halo of the Milky Wa y and 
to the halos of z ~ 3 LBGs l)Adelberger et al.ll2005l ). this 
should not be interpreted as supporting that LBGs are the 
precursors of galaxies like the Milky Way or z ~ 3 popula- 
tions. Linking a population of galaxies at one redshift with 
another popul ation at another is fraught u ncertainty (see 
for example, iMoustakas fc Somerv illc 2002). If these galax- 
ies are the progenitors of galaxies like the Milky Way, this 
would require an evolutionary path whereby the cluster- 
ing amplitude and halo mass does not grow, the bias de- 
creases, but the baryonic mass would have to increase by 
over an order of magnitude. I n a galaxy conserving model 
for the growth of structure IjFrvl 1 19961 ) and if we think 
of the estimated bias of the LBG population as z~5 as 
its bias at birth, would suggest that bias would decrease 
with dec reasing redshift while the correlation sca le would 
increase (|Frvlll996l : IMoustakas fc Somervillell2002h . In such 
a model, the evolution of both the bias and correlation 
scale suggest that the z ~ 5 LBGs could be the progeni- 
tors of the red galaxies at z ~ 1 and massive early type 
galaxies at z«0 jCoil et al] 12004 IZehavi et al] 120051 ). The 
con clusions would be s i milar in hierarchical merging m od- 
els (jOuchi et al ] |2004bl : IMoustakas fc Somervilld l2002h . In 
either model, it is difficult to say that the z ~ 3 and z ~ 5 
are directly related (except for the large uncertainties in the 
clustering and bias estimates for the z ~ 5 LBGs). 

In the hierarchical merging scenario, phase-space den- 
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sity arguments suggest that massive, present-day elliptical 
galaxies have formed from the merging of high redshift disk 
galaxies, w hich are smaller a nd denser than disk galaxies 
seen today l|Mao fc Molll99ct ). If the z~5 LBGs represent 
the early stages of the formation of present day bulges and 
spheroids they should have core phase-space densities that 
are consistent with their present-day descendants. A signifi- 
cantly lower phase space dens i ty wo uld r ule out a pro g enitor 
scenario. Following ICarlberd (| 19861 ) and iMao fc Mol (|l998T ) 
the core phase-space density of the typical z ~ 5 LBG in our 
sample is ~10~ 6 M© (pc km s -1 ) -3 . This is consistent with 
the core phase-space densities of massive, bright pr esent day 
ellipticals and the central region of the Milky Way (|Carlber el 
1 19861 ; IWvselll998l ; IMao fc Mo|[l998l ; lAvila-Reese et al.ll2005h . 
While we cannot constrain the evolutionary path of these 
LBGs, this similarity between the core phase-space densi- 
ties of our z ~ 5 LBGs and present day bulges and ellipti- 
cals, strongly suggests that the central population in these 
systems is already in place at high redshift. 

The young z ~ 5 LBGs may therefore represent the for- 
mation of the central regions of present-day massive galaxies 
at the earliest stages of their evolution. If so, their abun- 
dance at a look back time of 12-13 Gyr is consistent with 
the ag e of the do minant stellar population in the Galactic 
Bulge (|Richl2005l ) and its pr oposed formation in short bursts 
1-2 Gyr after the Big Bang (jZoccali et al.ll2006l L 



6.2 Enrichment of the intra- and inter-galactic 
medium 

6.2.1 Winds 

These galaxies are young relative to the age of the Universe 
at the epoch they are observed (~1.2 Gyr at the median red- 
shift). The young ages derived for the majority of z ~ 5 LBGs 
imply that these systems began appreciable star formation 
at Zf ~ 6 with the bulk of their stars formed at z < 6. As 
such, if reionisation w as complete at z~ 5.8 (jBecker et alj 
l200ll ; iFan et~al]l200lh . the intense star formation episode 
studied here could only significantly contribute to the end 
of reionisation and these galaxies are unlikely to be the per- 
petrators. However, these young galaxies do have an impact 
on the high redshift IGM. These compact galaxies are ex- 
periencing widespread and recent star formation at prodi- 
gious rates. The UV emitting regions are small (~1 kpc) 
and consequently these z ~ 5 LBGs have high rest-frame UV 
surface brightnesses. Given the typical star formation rates 
inferred from the SED modelling, their star formation rate 
surface densities (star formation rate per unit area) are sev- 
eral tens to hundreds Mq yr _1 kpc -2 . This is far higher than 
observed star formation rate surface densities of local galax- 
ies that are known to be driving vigoro us gaseous outflows 
(>0.1M Q yr" 1 kpc" 2 ; iHeckmanl 12001! ) . Thus, young z ~ 5 
galaxies host strong outflows, expelling the metals created 
by stellar nucleosynthesis into the surrounding medium. As 
such, these vigorous starbursts must contribute to the early 
chemical enrichment of the IGM which is supported by the 
detection of ions of m etals such as C arbon in the IGM at red- 
shifts as high as ~ 5 jSongaila 2005, and references therein). 



6.2.2 Galaxy-scale metal-mixing & metal-free star 
formation 

The first stars in the Universe must have been born from 
gas of primordial abundance (so-called population III stars) . 
Recent detections of the signatures of population III stars 
in the composite spectrum of z ~ 3 - 4 Lyman break galaxies 
suggest that such a population comprises 10-30% of the stel- 
lar mass, implying that galaxy- wide mix ing of metals is inef- 
ficien t on time-scales of a billion years (|jimenez fc Haimanl 
2006fl Most models of galaxy formation have assumed ef- 
ficient intra-galactic mixing of metals and do not predict 
metal -free star formation at redshifts significantly below 
z~5 jYoshida et alj|2004 ). 

While the star-formation intensities suggest that our 
young galaxies are likely to drive outflows and pollute the 
IGM with metals, and presumably their own interstellar me- 
dia, their very young ages imply that this material does not 
reach far beyond galactic scales. We can quantify this state- 
ment by making some simple assumptions about the out- 
flows driven by these high intensity star-forming galaxies in 
analogy with local starburst galaxies. 

The temperature of the outflowing gas can be estimated 
by simple energetic arguments (i.e., energy injected rate of 
the intense star formation equals the rate at whic h material 
- both ambient and stell ar ejecta - is thermalised; IHeckmanl 
l200ll ; lMoran et aljfl999h 

T wind » 0Am„-^- = 5.4 x W^f- 1 K (2) 
kM M10 

where mn is the mass of hydrogen, E is the energy in- 
jection rate, and M is the mass outflow rate, k is Boltz- 
mann's constant, and fi oa ding is the fraction of hot gas 
that results from "mass loading" due to sweeping up and 
entrai ning ambient interste llar material. Using the mod- 
els of iLeitherer et all (|l999T ). at the median age and star- 
formation rates we have estimated form fits to the SEDs 
of the robust sample of LBGs, the typical energy and mass 
injection rates are about 10 43 ergs s _1 (=£43) and about 
10 Mq yr -1 (=Mio), respectively to which we have scaled 
equation 16.2.21 

Similarly, equating the kinetic energy of the gas to the 
energy injection rate yields an estimate of the outflow ve- 
locity: 

W = (2§ )^ 2 ~ 1800/-/L 9 (# 1 ) 1/2 kms- 1 (3) 
M M10 

The mass loading of local starburst winds has been es- 
timat ed to be about a factor of a few to 10 l|Moran et al.l 
1999). This would decrease the wind velocity accordingly. 
If we assume that winds are loaded by about a factor of 5, 
then the wind velocity is «800 km s _1 . In the typical age 
of the LBGs, the wind would traverse ~15 kpc. Therefore, 
in principal, it is possible for the wind to span and enrich 
the ambient interstellar medium with metals. However, we 
have assumed that the mass loading is a simple drain on 

8 Note: our assumed 0.2Zq correspon ds very well to the combine d 
Z~ OZq and Z~O.4Z0 mixed model of lJimenez fc Haimanl ||2006|) . 
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the energy of the wind material itself. In reality, the am- 
bient ISM will be an inhomogeneous medium with clumps 
of higher density material. The diffuse ISM will get swept 
up while the clouds will be entrained and accelerated by the 
ram pressure of the wind material itself. In such an situation, 
it is important to consider both how a continuous medium 
and clouds are accelerated in the outflowing gas. Following, 
e.g.. I Dyson fc Williams! (|l980T l. expansion speed of the shell 
and energy content of the bubble are related as: 

VsHM ~ 240S43 75 n" 1/5 t~ 2/5 km s" 1 , (4) 

where E43 is the (constant) energy injection rate in 
units of 10 43 ergs s~\ no — 1.0 cm~ 3 is the ambient ISM den- 
sity in cm" 3 and tr is the injection time in units of 10 Myr. 
At such velocities, swept-up ambient ISM will traverse about 
5 kpc over the age of the typical LBG in our sample. Clouds 
accelerated in the wind will be dispersed over a similar size 
scale. Since this is only a factor of a few larger than the half 
light radii of the galaxies in our sample, given the strong col- 
limation of superwinds observed in the local Uni verse (e.g., 
iHeckman et~alj Il990l ; iLehnert fc Heckman|[l99r} ). it is un- 
likely that the metals in the outflow will be mixed signifi- 
cantly with the infalling or ambient gas. 

While even the swept-up ambient ISM is unlikely to 
be efficiently mixed by the pressure and outflow of the in- 
tense star-formation, another interesting time scale to in- 
vestigate is the cooling time of the hot superwind gas. If 
the gas were to have a long cooling time, then it would 
be unavailable for star-formation. This would effectively se- 
quester the metal enriched gas from the stellar mass loss 
from being incorporated into star-formation over a cooling 
time scale. We estimated previously that the wind mate- 
rial itself, modulo the mass loading, will have a temperature 
of « 5 x 10 7 K. This plasma will also be very rare. We can 
crud ely estimate the den s ity of the gas using the relationship 
from I Chevalier fc Cleed l|l985h . namely, 

n e>wind w M^E-^R- 2 (5) 
« 6 x 10 _2 M 1 1 6 s S43°- 6 ^ 2 cm" 3 (6) 

where we have scaled the radius to lkpc, Rikpc, ap- 
proximately the half-light radius of the typical LBG in our 
sample and the area over which energy is being injected 
by massive stars. If the wind that has reached large dis- 
tances from the LBG is mass loaded, the density will be of 
course higher (likely be a factor of a few to 30). The various 
ranges of density, depending on radius and mass loading is 
of-order 10 -2 to 10 -4 , similar to what has been estimated 
in local starburst galaxies. Using this crude estimate, we 
can investigate the likely cooling time of the wind using 
standard cooling arguments. For an equilibrium plasma at a 
temperature of 10 6 7 to 10 7 T K and density of 10~ 2 tol0" 4 , 
the cooling time is, roughly 10 s years to a few x 10 9 years. 
Over that duration, the star-formation will be able to build 
up about 10-30% of the mass of LBGs at z ~ 3 and do 
so with much of the metals generated by the intense star- 
formation at higher redshifts locked in hot gas that takes 
about 1 Gyr to cool. Of course, the velocities estimated are 
also to likely be approximately or greater than the escape 
velocities of the dark matter halos and much of the enriched 



material may escape. These tight timing constraints are con- 
sistent with evidence of both a significant fraction of popu- 
lation III stars in lower redshift Lyman break galaxies and 
inefficient metal-mixing o n intra- and inter-galactic scales 
(IJimenez fc HaimanllioOrj '). 



7 SUMMARY 

We present the properties of a robust sample of LBGs at 
z ~ 5 selected to %ab < 26.3 in the Chandra Deep Field 
South. On average, LBGs at z ~ 5 are ~10 times less massive 
(~10 9 M Q ) and are significantly younger (< 100 Myr) than 
similarly luminous LBGs at z ~ 3. While LBGs with such low 
masses and young ages are not unknown at z ~ 3 they are far 
less common, such systems comprise >70% of LBGs at z ~ 5 
and in contrast only 30% at z ~ 3. Their short duty cycles 
suggest that the z ~ 5 population must be highly stochastic 
and that samples of z ~ 5 LBGs may be highly incomplete, 
with only ~l-in-13 LBGs being detected. Considering this 
implies that up to ~1.5% of the present-day stellar mass 
density can be accounted for by luminous LBGs (c.f. 15% 
a Gyr later iRudnick et al.ll200q) . The abundant fraction of 
young and moderately massive z ~ 5 LBGs are likely to be 
experiencing their first (few) generations of large-scale star 
formation and are accumulating their first significant stellar 
mass. Their formation redshifts (z~6-to-7) suggests these 
galaxies could have contributed only to the end of reionisa- 
tion. They are extremely compact as measured from the UV 
ACS images (ri/j ~ 1 kpc) giving rise to high stellar mass 
and star formation rate surface densities. The stellar mass 
surface densities and core phase-space densities are compa- 
rable to the spheroidal components of present day L* galax- 
ies and suggests these systems are the progenitors of early- 
type galaxies or bulges, favouring the inside-out galaxy for- 
mation scenario. The high star formation rate surface den- 
sities implies these sources are driving winds that enrich the 
inter and intra-galactic media with metals. Their young ages 
are consistent with inefficient metal-mixing on galaxy-wide 
scales. Therefore it is plausible that these galaxies contain a 
significant fraction of metal-free stars as has been proposed 
for z ~ 3 LBGs. 
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APPENDIX A: ALTERNATIVE MODELLING 
ASSUMPTIONS 

To assess the sensitivity of our characteristic ensemble prop- 
erties to the model assumptions, we fit the SEDs of the z ~ 5 
LBGs with different sets of input parameters, varying the 
star formation history (SFH), extinction law and metallic- 
ity. As we are interested in the ensemble properties, we fo- 
cus on the effects of variations on the median properties of 
the robust sample of z ~ 5 LBGs. Degeneracies between the 
model parameters make it difficult to discriminate between 
SFHs, extinction laws or metallicities based solely on the re- 
duced chi-squared values (x„). One must also keep in mind 
that there are uncertainties in the stellar evolution, spectral 
libraries, and dust properties themselves, especially impor- 
tant at sub-Magellanic m etallicities (see the discussion by 
iBruzual fc Charlotl l2003) as well as regarding their validity 
at z ~ 5. Differences in goodness-of-fit, even systematic, may 
reflect in part template mismatches. 

We discuss the effects of varying the assumed parame- 
ters in the following subsections and Table |A"T1 summarises 
some of the key variations. However, in view of all the un- 
certainties involved, the numbers quoted in the following 
should not be taken at face value but rather as indications 
of overall trends. The interloper fraction (low z galaxies as 
well as star/QSO candidates) for the initial sample of f09 
candidate z ~ 5 LBGs is roughly 40% in all cases; the me- 
dian z p h for the robust sample is always around 4.8 except 
for dust-free models with star formation time-scales longer 
than ~ 50 Myr, for which it increases by about 0.2. Only 
models without dust extinction could affect our conclusion 
of young ages and could make the ensemble as old or older 
than ~ fOO Myr. However, such models are also those that 
lead to the poorest fits compared to any other set of in- 
put parameters, although they cannot be ruled out from a 
purely statistical \ 2 argument. Typical stellar masses for the 
ensemble of a few x f 9 Mq are derived for all but the lowest 
metallicities explored, Z = 0.002 - 0.02 Z . 

Al Star Formation History 

In addition to our adopted constant star formation rate 
(CSF) scenario, we considered a suite of exponentially de- 
clining star formation rates (SFRs) with e-folding time- 
scales (r) of 10, 30, 50, 100, 300, 500 Myr and I Gyr, as 
well as the case of a single stellar population (SSP) formed 
in an instantaneous burst. 

Adopting the 0.2 Z models and the SMC extinction 
law, the median best fit age is a few tens of Myr, for all 
SFHs except for the SSP where it decreases by a factor of 
~2. This behaviour is not surprising given the degeneracies 
between age and SFH: the best-fit age gives a measure of 
the time elapsed since the bulk of stars that dominate the 
observed SED was formed, so the derived ages tend to be 
younger for models with shorter star formation time-scales. 
Since most of the LBGs in the robust sample have ages of 
a few tens of Myr, any SFH with time-scale of that order 
or longer is effectively equivalent to a CSF model and the 
derived ages are similar. The best-fit extinction changes lit- 
tle, with median Av values of a few tenths of magnitude 
for t > 30 Myr down to mag for shorter time-scales. The 
median SFR decreases by a very modest ~ 10% between 



CSF and r = 30 Myr models, and drop by a factor of 2 for 
t — 10 Myr (by definition, there is no star formation for 
a SSP with non-zero age). The median stellar mass shows 
no systematic variations and remains roughly constant with 
star formation time-scale, the largest difference amounting 
to ps 17%, smaller than the dispersion among individual ob- 
jects in the sample. 

The trends of lower median best-fit ages, Av values, and 
SFRs with star formation time-scales are qualitatively the 
same irrespective of the metallicity or extinction law. The 
exception is the median age for dust-free 0.002 Zq models, 
which increases from 3 Myr for all SFHs to 6 Myr for an SSP; 
this difference is however at the limit of the age sampling in 
our modelling and is not very significant. The median SFRs 
drop by « 15% up to a factor of sa 5 from CSF to r = 10 Myr 
models, depending on the metallicity and extinction law. 
For the stellar masses, the median values vary by a factor of 
1.5 or less between different SFHs for any given metallicity 
and extinction law considered, and up to a factor of 1.8 for 
models without dust extinction. 

A2 Extinction Law 

Our choice of an SMC- type extinction law was motivated by 
consistency with the adopted metallici ty of Z = 0.2 Z . As- 
suming the starburst extinction law of ICalzetti et al.l (|2000l . 
hereafter Calzetti law) instead, which is widely applied in 
high redshift studies, the median best-fit Av values increase 
by a large amount up to about 1.2 mag (0.6 mag for an SSP) 
and at the same time, extremely young median ages of only 
a few Myr are derived. The Calzetti law is much greyer (i.e., 
shallower) in the rest-frame UV than the SMC law, so that 
higher levels of extinction are allowed but this then requires 
younger ages to fit the very blue SEDs of our LBG sample. 
As a consequence of the younger ages and higher Av values, 
both the median stellar mass and SFR increase, by factors 
of 2 — 3 and 30 — 60, respectively. 

The properties and geometry of interstellar dust at 
z ~ 5 is currently an open issue, and the application of ex- 
tinction laws derived for low redshift galaxies to high redshift 
systems is questionable. Low to moderate extinction would 
be expected for young systems in which it would be unlikely 
that significant amounts of widespread dust have had suf- 
ficient time to form. Moreover, the detection of Lya in the 
spectrum of some of these systems argues against strong at- 
tenuation by dust. There is empirical evidence t hat dust in 
z ~ 2 — 4 LBGs exhibits SMC-like characteristics (Vijh et al.l 
l2003h . At higher z and given the young ages inferred for our 
LBG sample, it may be that dust from Type II supernovae 
dominates the extinc tion, as proposed for z~6 quasars by 
iMaiolino et al.l l|2004h . However, in the absence of determi- 
nations valid for star formation dominated z ~ 5 LBGs, 
the SMC extinction law is a reasonable assumption. It also 
leads to modest extinction, older ages, and lower SFRs that 
are more plausible than the extremely young ages of a few 
Myr, and the high Av ~ 1 mag and SFRs ~ 1000 M yr" 1 
obtained with the Calzetti extinction law. 

Models where we allowed extinction to be a free pa- 
rameter fit the data better than models without attenua- 
tion by dust, even though the required best-fit extinction is 
modest for our adopted metallicity and SMC extinction law 
(Av ~ 0.3 mag). This is reflected in the overall systemati- 
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cally higher Xn values for the zero-extinction case. Strictly 
speaking, the differences are not statistically significant so 
that models without dust extinction cannot be conclusively 
excluded. However, given that there is spectr oscopic evi- 
dence that the metallicity of LBGs is 0.2 Zq (|Ando et al.l 
120041 ). a dust-free solution is highly unlikely. 



A3 Metallicity 

We explored the effects of metallicity with models for Z — 
0.002, 0.02, 0.2, 0.4, 1 Z Q (i.e. those available for the Bruzual 
& Chariot models based on the Padova 1994 evolutionary 
tracks and BaSeL 3.1 spectral library). For this comparison, 
we chose the Calzetti extinction law for all metallicities. For 
any of the SFHs we considered, the best-fit ages are slightly 
older at higher metallicities. The highest median extinction 
(Av ~ 1 mag) is derived for Z = 0.2 and 0.4 Zq models. The 
median stellar masses decrease by a factor of 4 — 5 between 
0.002 Zq and Zq while the median SFRs vary by a factor of 
3 — 8 and are typically lower at both ends of the metallic- 
ity range explored. The overall trend for the median stellar 
masses is preserved for models without dust extinction, but 
those for the ages and SFRs are reversed. 

The typical stellar masses of our robust sample LBGs 
are thus most affected by changes in the metallicity, but 
by only a factor of 5 or less. The median masses for our 
adopte d 0.2 Zq metallicity , supported by spectroscopic re- 
sults of lAndo et al.l (|2004h . lies at the middle of the range 
inferred between 0.002 Z Q and Zq for any given SFH. 



APPENDIX B: LIMITS ON THE 
CONTRIBUTION OF AN OLD POPULATION 

The SEDs of all our z ~ 5 LBGs are very blue in the rest- 
frame UV, indicating the presence of on-going or very recent 
star formation. The best-fit stellar masses (like all other de- 
rived properties) represent those of the young and luminous 
populations that dominate the integrated rest-frame UV to 
NIR light. By selection, we would not expect to have picked 
up old passively evolving galaxies. Therefore it is unlikely 
that these luminous LBGs contain a significant, older, less 
luminous component that significantly contributes to the 
stellar mass. Moreover, the z ~ 5 LBGs in our sample are 
extremely compact in the GOODS-ACS data (ri/j ~lkpc), 
and given that they are selected at an i-band magnitude 
significantly brighter than the depth of these images, any 
old stellar population would have to follow an extremely 
broad distribution to remain undetected. If any significant 
old population is distributed over the same extent as the UV 
luminous galaxy, then this would imply an extremely high 
mass surface density (> 10 10 Mq kpc~ 2 ), in place already at 
z ~ 5. Such high mass surface densities are known only in the 
inner regions (<kpc-scale) of massive cD galaxies, and they 
are not for typical starburst galaxies or z ~ 3 LBGs. This 
unlikely high surface density strongly argues against a sig- 
nificant fraction of old stars contributing to the galaxy mass, 
and suggests that the stellar mass derived from the stellar 
synthesis modelling is representative of the total mass of the 
galaxy. 

Nevertheless, we can place extreme upper limits on the 



mass contribution of such an old underlying population fol- 
lowing a similar procedure to that used for LB Gs at z ~ 3 
l|Papovich et al.ll200ll ; IShaplev et aDl200lL l2005h . This con- 
sistent approach allows us to test whether such an old contri- 
bution can explain the mass difference we have found which 
has been discussed in Section [5] For this we consider an ex- 
treme 'worst case' scenario, requiring the coincidence of a 
SSP which, compared to scenarios of extended star forma- 
tion, maximizes the mass-to-light ratio of a "hidden" compo- 
nent for a given age. For this we assume the true photometry 
of the source is in fact at the brightest la limit of our photo- 
metric uncertainties. These uncertainties already include a 
generous 10% to account for flux calibration errors. There- 
fore, the limit that results is extreme and while there is a 
chance that it can apply to an individual source, it is highly 
unlikely to be valid for the population as a whole. 

We consider the case of a passively evolving component 
with an age of 500 Myr, the median for the one-third of 
our robust sample older than 100 Myr, and comparable to 
the oldest a ges for z ~ 5 LBGs reported by ot her authors 
as well (e.g. lEvles et al.ll2005l : IVan et alJl2006bl ). The SED 
of a 500 Myr old SSP peaks in the rest-frame optical, in 
particular just red- wards of the Balmer/4000 A break, and 
so would contribute most to the observed fluxes in the K s , 
3.6 /im, and 4.5 /im bands. We thus determined the upper 
limit for an old component by adding to the nominal best- 
fit SED for each object the maximum contribution from a 
500 Myr old SSP allowed by the lcr uncertainties on the 
K s , 3.6 /xm, and 4.5 /im photometric points. The typical 
maximum contribution corresponds to an increase in flux 
by a factor of 2.5 and translates into a factor of m 7 times 
more stellar mass0. 

However this factor must not be taken at face value. 
The assumption of a 500 Myr - to maximally-old passively 
evolving population formed in an instantaneous burst, made 
to maximise the mass-to-light ratio, is more extreme than 
any detected population. Any scenario with an old compo- 
nent that has been forming stars for some non-negligible 
period of time, even if short, will reduce this upper limit. 
And considering our conservative photometric uncertain- 
ties, an upper limit to the mass contributions for our z ~ 5 
LBGs will be considerably less than thi s. Similar factors have 
been determined fo r z~3 LBGs (3-5 IPapovich et al.l l200ll : 
IShaplev et afll200ll L and therefore such a potential contri- 
bution from old stars would be insufficient to make the mass 
ranges shown in Section Oof the sample overlap. 
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Table Al. The effects on the derived properties of SFR, mass, age and extinction of our z~5 LBGs with a selection of varying input 
parameters are summarised in this table. 
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